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Abstract - We exactly calculate the reduced density matrix of matrix product states (MPS). Our 
compact result enables one to perform analytic studies of entanglement in MPS. In particular, 
we consider the MPS ground states of two anisotropic spin chains. One is a g-deformed Afneck- 
Kennedy-Lieb-Tasaki (AKLT) model and the other is a general spin-1 quantum antiferromagnet 
with nearest-neighbor interactions. Our analysis shows how anisotropy affects entanglement on 
different continuous parameter manifolds. We also construct an effective boundary spin model 
that describes a block of spins in the ground state of the q-deformed AKLT Hamiltonian. The 
temperature of this effective model is given in terms of the deformation parameter q. 
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Introduction. — We study entanglement in the 
ground state of spin models. Specifically, we consider spin- 
1 chains with nearest-neighbor interactions. We begin by 
analyzing the reduced density matrix pi of a block of se- 
quential spins within the chain. The reduced density ma- 
trix is fundamental in the study of entanglement. It is 
used to define the Renyi and von Neumann entanglement 
entropies 

logtrp? 



Sr(o) 



l-a ' 
-tr(p^log^) 



lim S-R(a) 



(1) 
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These quantities are frequently used measures of entan- 
glement []}j4]- Also, the eigenvalues of pe are used to 
construct the entanglement spectrum of the block [5j|8j. 
This spectrum describes the mixed state of the block as 
if it were governed by an effective Hamiltonian at some 
temperature. 

In this paper, the reduced density matrix of translation- 
ally invariant matrix product states (MPS) [9 18 is calcu- 
lated exactly. This class of MPS includes one-dimensional 
valence-bond-solid (VBS) stat es [l2}|l6] and ground states 
of anisotropic spin chains [17||18| . In particular, we con- 
sider the MPS ground states |MPS) of two nonequivalent 
anisotropic generalizations of the Affleck-Kennedy-Lieb- 
Tasaki (AKLT) model (L2l[l^. Both are spin-1 chains 



described by Hamiltonians of the form H 
The local Hamiltonian h 



j.j+i acts on the Hilbert space 
of neighboring spins j and j + 1. Additionally, it anni- 
hilates the frustration-free ground states of these models, 
/i iJ+ i|MPS) = 0. 

The first generalization we consider is a g-deformed 
AKLT (AKLT g ) model that is SU g (2) invariant fl7l[l9 



The second model is a general Hamiltonian with U(l) 
symmetry introduced by Klumper, Schadschneider, and 
Zittartz (KSZ model) [18). The KSZ Hamiltonian is in- 
variant under (i) lattice translations and reflections, (ii) 
spin rotations about the longitudinal z-axis, and (iii) spin 
reflections about the transverse plane Sf — > —S*. We 
discuss each in more detail in the following subsections. 

Entanglement in the VBS ground state of isotropic 
AKLT models has been studied in the literature 20 25 



However, the question of how anisotropy modifies entan- 



glement in these states has received less attention 26 ■ 28 



We address this matter by calculating the entanglement 
spectra and entropies in AKLT 9 and KSZ ground states 
as functions of anisotropy parameters. 

Reduced density matrix of MPS. We start with 
a pure matrix product state |MPS). The density matrix 
of the whole state is p = |MPS)(MPS|/(MPS|MPS). We 
then partition the system into a block of t sequential spins 
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and its environment E. The reduced density matrix pe is 
the partial trace of p over the environment, pi = tr^ p. 

The MPS representation of a periodic chain of L iden- 
tical spins is 



IMPS) =tr[ gl .g 2 ..... gi ]. 



(3) 



The gj- are D x D matrices (D = 2 for the examples con- 
sidered below). The trace here is taken over the auxiliary 
matrix space (not E). The elements of g- are 



(g,W 



(m)|m) 



(4) 



The set {|m)j} is a complete orthonormal basis for the 
Hilbert space of the spin at site j and the coefficients 
A a /j(m) are independent of the site index. Due to trans- 
lational invariance, we drop the site label j whenever pos- 
sible. We denote the matrix dual to g as g with elements 
(g)a/3 — Em ^as( m )( m l- Here, the coefficients are re- 
placed by their complex conjugates and the kets are re- 
placed by the corresponding bras. In Eq. (|3| the matrix 
multiplication (•) involves tensor products of vector matrix 
elements, i.e., 



(gj-gj+iW = A a p(m)A^(n) \m)j ® \n) 



3+1- 



(5) 



(3mn 



The dual (gj-Ej+ijaj is defined analogously. For prod- 
ucts of g matrices denoting a block of sequential spins we 
introduce an abbreviation 



(gj-gj+i"---"gfW = \aa';j,j' 



(6) 



Thus, it is sufficient to identify the states of boundary 
spins to specify the state of a block. 

Let us construct a transfer matrix G = g eg g that is 
useful for calculating state overlaps (scalar products) and 
correlation functions: 

(G)a 7 ,/3S = (g)af)(g)-,6 = ^ Kfli™) ^M m )- (7) 



For example, the square of the norm of |MPS) is 
(MPS|MPS) = J2 aa ,(aa';l,L\aa';l,L) = trG L . 
The density matrix is therefore 



•gL]tr[gi'----gL] 



(8) 



\aa'; 1, £) (G L ~ e ) a ,^, a/3 ((3(3'; 1, l\ ^ ^ 



tvG L 

while the reduced density matrix is 

-L 

\aa ; ±, n) \\ 

Pi 



aa'13/3' 



tiG 1 



We chose the block to extend from site 1 to site £ with- 
out loss of generality because of translational invariance. 
In this form, pi clearly acts nontrivially only in the sub- 
space spanned by the block state vectors {|aa'; The 
dimension of this subspace is at most D 2 . 



Let us consider transfer matrices with the symmetries 
(G) a p, a >p> = {G) a 'p>,ap = {G)/3a,p'a>- This is a weak 
requirement because it means that the scalar product 
(aa'; 1, n|/3/3'; 1, n) is invariant under lattice reflections 
(e.g., 'flipping' the ring over). We construct a symmet- 
ric overlap matrix K(n) that is related to the n th power 
of G by (K(n)) aa , )/9i g, = (G n ) a p, a >p>. With this definition 
we write pe as 



Pc 



aa'13/3' 



trG 1 



The indices are now matched so that we can express pi as 
a product of matrices. Suitable similarity transformations 
within the space spanned by {|aa'; 1,^)} finally gives 



K(L-£)K{£) 
tr[K{L - £)K(£)Y 



(11) 



This formula has some notable features. First, it is a 
general expression that is valid for a large class of MPS. 
Also, we find that pi has a small number of nonzero eigen- 
values, rankpf < D 2 [29j[30]. Furthermore, it is evident 
that pg and pl-i are isospectral. Thus, the entanglement 
entropies of the block and environment are the same. 

AKLT 9 model. — The AKLT g Hamiltonian is given 
in the Appendix. It is hermitian for real q and has the 
symmetry of the SU g (2) quantum group 31 32 . This 
deformation of SU(2) symmetry appears naturally in the 
anisotropic spin- 1/2 XXZ Heisenberg model [33] . 

The generators of the SU g (2) quantum algebra are 



Ej &j' ■ They satisfy the commutation relations 



q 2S> _ q -2$» 



q-q 



l - r - , |1S* I S ± ]=±S ± . (12) 



The q-deformed spin-1 operators at site j are S| = and 



(13) 

Here 5?' are undeformed spin-1 operators at site j. 



Eq. ( 13 1 has similarities with Jordan- Wigner transforma- 
tions. 

The entanglement spectrum does not depend on the sign 
of q and is invariant under the transformation q — > q . 
Thus, we will only consider q € (0, 1]. The isotropic AKLT 
model is recovered at q = 1. 

The unique ground state of the periodic AKLT 9 model 
is a VBS 9 state. Its MPS representation was constructed 
in Ref. [17]. The properties of this state and its spin- 
spin correlation functions have been studied 17 34 . This 



approach was later extended to obtain the spin-spin cor- 
relators of higher integer spin-S* AKLT g chains 



3(> 



The local Hamiltonian hjj+i is a projector onto the 
subspace of the g-deformed spin-2 quintuplet formed by 
adjacent spins [17]. The MPS form of the VBS g ground 
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Fig. 1: (Color online) (a) In the double scaling limit the re- 
duced density matrix of |VBS 9 } has four degenerate eigenval- 
ues at the isotropic q = 1 point. The two middle eigenvalues 
are always degenerate (bold line). In the classical limit q — > 
the ground state is a product state and the only nonzero eigen- 
value is unity, (b) The Renyi entropy Sr vanishes as q — > 0. It 
saturates to the maximally entangled value Sr = log 4 at the 
AKLT point q = 1. The von Neumann entropy (bold line) is 
obtained at lim a ^i Sr. 



state is obtained by requiring hjj+i to annihilate the ele- 
ments of the matrix g ? - Sj+i- This condition leads to 



g 



Jq~+q~ 1 \-) 



V'? + ^ 1 l+> 

-q\0) 



(14) 



where the vector elements are eigenstates of S z . 
yields the MPS representation 

|VBS 9 )=tr[ gl .g 2 .....g L ]. 



(15) 



Since g is a 2 x 2 matrix, the reduced density matrix has 
at most four nonzero eigenvalues. 

In the limit q — > 0, the AKLT^ Hamiltonian is domi- 
nated by (classical) Ising-type interactions. In this case 
the g matrix has only one diagonal element |0). The re- 
sulting ground state is a product state • \0)j describing 
a magnet polarized in the transverse direction. In this 



classical limit all spins are in the Si- = state. Hence 
any block in the chain has zero entropy. 

The transfer matrix ^ for the |VBS g ) state is 



/ 



g®g = 



q + q- 
-10 
0-1 



\q + q~ 







The eigenvalues of this matrix are {A, - 



dominant eigenvalue is A 
matrix K(£) is therefore 



l + q z + q 



q 

-l. -l 
2 >3. 



(16) 



—1} and the 
The overlap 
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q + q' 
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- diag{q \ 1, l,q} 
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q + q- l \ 



-1 

q- 1 J 



(17) 



Eq. ( 11 ) gives the reduced density matrix 



Pi = 



K(L-£)K(t) 
A L + 3(-l) L ' 



(18) 



In the double scaling limit of an infinite block £ — > oo 
in an infinite chain [L — £) —¥ oo, the reduced density 
matrix becomes diagonal. Thus, the block states \aa'; 1, 1) 
are orthogonal to each other. The eigenvalues {pi} of the 
reduced density matrix pe~^oo are 



7 T2 



1 



Pl,i 



(q + q- 1 ) 



P2 = P3 



(q + q- 1 ) 



M2' 



-►oo. (19) 



We discover an important consequence of g-deformation: 
The degeneracy of the entanglement spectrum changes be- 
tween the classical and isotropic AKLT points. This result 
is depicted in Fig. [TJa). 

To obtain an intuitive picture for the mixed state of 
the block we write pe^oo — e^ l3H " / tr e~ l3H " . H e is an 
effective Hamiltonian and 1//3 an effective temperature. 
The eigenvalues of H e form the entanglement spectrum [5] . 
Doing so gives the effective temperature 1//3 = 1 / | In ^ | 
and effective Hamiltonian 



Hp 



a 1 a e 



(20) 



This where are Pauli operators at site j. Thus, pe^co 
describes a thermal ensemble of two spin-l/2's at the 
block boundaries with Heisenberg (XX) interaction. The 
anisotropy parameter q determines the effective boundary 
temperature T e = 1/ \ lnq\. For the original AKLT model 
(q = 1) the effective boundary spins are in a maximally 
mixed state (T e — > oo), while in the classical limit q — » 
they are in a pure state (T e = 0). 

This result for the effective Hamiltonian is consistent 
with the area law for gapped models 37|38 . It is similar to 
the effective boundary spin chain proposed for 2D AKLT 
dels 0139] 



However, in the AKLT 9 chain the effective 
boundary spin interaction is long-ranged and exists for 
arbitrarily long blocks. The long range of this interaction 
follows from the non-local nature of the SU ? (2) symmetry 



( 13 1 of the model. 



The Renyi and von Neumann entropies in the double 
scaling limit are 



2 q* 
log- 



q 



~ l-a iU& (q + q- 1 ) a ' 



Si^^logiq + q-'r- 



q + q- 



(21) 
(22) 



We plot these entropies in Fig. [TJb) and observe the 
effects of g-deformation even for infinitely long blocks. 
Anisotropy reduces entanglement entropy from its max- 
imum value Sb. = log 4. This maximum is reached at the 
AKLT point (^). 

For blocks of finite length I < oo in an infinite chain L — > 
oo, the overlap matrix K(^) has off-diagonal terms. In this 
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- h<p.4) 
■ log(3) 




(a) 



(b) 



Fig. 2: (Color online) The von Neumann entropy decreases 
away from the isotropic AKLT points q = 1 for |VBS 9 ) (a) 
and a — 2 for |KSZ a ) (b). At the isotropic point finite-size 
corrections are largest for |VBS,j} and smallest for |KSZ a ). 



case, the eigenvalues of pt acquire finite-size corrections: 



q 2 + q 2 



PlA 



P2=PZ = 



2(-A)" 



2(q + q-l)2 

1 - (-A)- £ 
(q + q- 1 ) 2 ' 



± 



l-(-A) 

(q + q- 1 ) 2 



-2t 



A = 1 + q 2 +q 2 . 



(23) 



These eigenvalues are exact. The corrections decay expo- 
nentially as expected for a gapped system. At q = 1 we 
recover the eigenvalues 



l + 3(-3)" 



l-(-3)" 



Pl= V , P3=P3=P4 = ^^ , (24) 

for the isotropic AKLT chain [20]. Figure [2] a) shows how 
g-deformation reduces the von Neumann entropy. For a 
block consisting of one spin {I — 1) one eigenvalue is iden- 
tically zero. Thus, the maximum single-site von Neumann 
entropy is log 3. This value is reached at the isotropic 
point q = 1 where there is a uniform mixture of three 
spin-1 states. 

For large but finite blocks 1 <C I < oo and q ^ 1, an 



expansion of the eigenvalues ( 23]) gives the leading order- 
corrections ±( — l)^(f? + q^ 1 ) t 



£ /£. The characteristic 



length of these corrections is £ = 1/ ln(l + q 2 + q 2 ). This 
quantity is equal to the correlation length of the spin-spin 
correlators of the VBS, state [l7}[34j[36] . 

KSZ model. — The local KSZ Hamiltonian is given 

by 

hjj+i = a A 2 + a 1 (A j Bj + BjAj) + a 2 B 2 + a 3 Aj 



+ aiB^l + Bj) + a 5 [(S|) 2 + 
with a transverse interaction term A 



QX QX 

qz qz 



a 6 , (25) 
I qV qy 
and constants 



longitudinal interaction term Bj 

a{. Requiring hjj+i to have nonnegative eigenvalues and 
annihilate an MPS ground state |KSZ a ) leads to a sub- 
manifold of Hamiltonians with restrictions on the con- 
The correlation functions and low-lying 

but 



stants a* 



18 



excitations of this model have been studied 18 40 



its entanglement spectrum has not yet been investigated. 



We obtain the MPS form of |KSZ a ) from the g matrix 

(26) 



( |0> 
W5|- 



-a|0) 



where a > is an anisotropy parameter and a = sgna3- 
The corresponding transfer matrix is 



(27) 



The unique ground state is |KSZ a ) = tr[g x - g 2 - . . . • g L ] 
(periodic boundary conditions). It reduces to the isotropic 
VBS state at a — 2 and a = 1. 

The eigenvalues of for an infinite chain L 
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oo are 
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PlA = j 
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1 



±2 



P2 



1 



(28) 



We observe that the entanglement spectrum is the same 
for a = ±1. In the double scaling limit the eigenvalues of 
Pe-¥oo are four- fold degenerate pl~*°° = \. The block is 
maximally entangled with Sr, = log 4. The entanglement 
spectrum therefore corresponds to a four-level system at 
infinite temperature. 

Let us now consider blocks of finite length. The von 
Neumann entropy is a maximum at the isotropic point 
a = 2. This property is depicted in Fig. [2]b). For a 
block of one spin {1 = 1) one eigenvalue of pi vanishes and 
the maximum entanglement entropy is log 3. In the limit 
a — > 0, the |KSZ a ) ground state approaches the transverse 
ferromagnet |0)j. This is a (classical) product state 
with no entanglement. In the opposite limit a — > oo the 
reduced density matrix represents a uniform mixture of 
two degenerate Neel ordered states. In this limit the von 
Neumann entropy approaches log 2. 



Finite-size corrections to the eigenvalues ( 28 ) are expo- 
nential in £. The characteristic lengths of these corrections 
are £|| = 1/ln |(1 + a)/(l - a)| and £± = l/ln(l + a). 
These quantities are equal to the longitudinal (£y) and 
transverse correlation lengths of the spin-spin corre- 
lation functions 18 : 



[sgn(l — a)] x e 



-mil 



(29) 



(l-o)"' 

(SfSf) = -a(a + l)[sgn(-a)] e xe-^, l>2. (30) 

Conclusions. — In this paper we studied entangle- 
ment in the MPS ground states of the anisotropic AKLT g 
and KSZ models. We derived a compact formula for the re- 
duced density matrix in terms of the MPS transfer matrix. 
With this result we calculated entanglement entropies on 
continuous parameter spaces connecting the isotropic VBS 
state with classical states. 
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In the AKLT g model we found that q-dcformation re- 
duces entanglement in the |VBS g ) state even in the dou- 
ble scaling limit. Our analysis enabled us to construct 
an effective Heisenberg model for the boundary spins of 
infinitely long blocks. This important result identifies the 
deformation parameter |ln q\ 1 as an effective temperature 
for these boundary spins. 

Anisotropy in the KSZ model also reduces entanglement 
in the |KSZ a ) state for blocks of finite length. But unlike 
the |VBSq) state, the boundary spins of infinitely long 
blocks are maximally entangled at fixed anisotropy < 
a < oo. Thus, the effective boundary spins in the double 
scaling limit are at infinite temperature. 
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Appendix. — The AKLT g Hamiltonian is 

n = b J2i cS i ' S i +1 + [ S J ' S J' +1 

j 

+ \{l-c)( q + q - 1 -2)S z J S z ]+1 
+ l(l + c)(q- q - 1 )(S* +1 -S*)] 2 
+ l c (l-c)( q + q -i-2) 2 (S*S* +1 ) 2 
+ \c{l + c){q- q - l ){ q + q- 1 - 2)S*S* +1 

+ i( c _3)[( c _l+i(l + c )2)^5« +i 
+ 2( c -|(l + c) 2 )((5; +1 ) 2 + (^) 2 )] 

+ ( C _l) + I C (g3_g-2)(^ +1 (31) 

with c = 1 + q 2 + q - 2 and b=[c{c- l)]" 1 [l9j[34]. 
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